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ABSTRACT 

We present the results of 16 S'wi/t-triggered Gamma-ray burst (GRB) follow-up observations taken 
with the Very Energetic Radiation Imaging Telescope Array System (VERITAS) telescope array from 
2007 January to 2009 June. The median energy threshold and response time of these observations 
were 260 GeV and 320 s, respectively. Observations had an average duration of 90 minutes. Each 
burst is analyzed independently in two modes: over the whole duration of the observations and again 
over a shorter timescale determined by the maximum VERITAS sensitivity to a burst with a 
time profile. This temporal model is characteristic of GRB afterglows with high-energy, long-lived 
emission that have been detected by the Large Area Telescope on board the Fermi satellite. No 
significant very high energy (VHE) gamma-ray emission was detected and upper limits above the 
VERITAS threshold energy are calculated. The VERITAS upper limits are corrected for gamma-ray 
extinction by the extragalactic background light and interpreted in the context of the keV emission 
detected by Swift. For some bursts the VHE emission must have less power than the keV emission, 
placing constraints on inverse Compton models of VHE emission. 
Subject headings: astroparticle physics - gamma-ray burst: general 



^ Fred Lawrence Whipple Observatory, Harvard-Smithsonian 
Center for Astrophysics, Amado, AZ 85645, USA 

^ Department of Physics and Astronomy, Barnard College, 
Columbia University, NY 10027, USA 

^ Department of Physics and Astronomy, University of Cali- 
fornia, Los Angeles, CA 90095, USA 

* Santa Cruz Institute for Particle Physics and Department of 
Physics, University of Cahfornia, Santa Cruz, CA 95064, USA 

° Department of Physics, Washington University, St. Louis, 
MO 63130, USA 

® School of Physics and Astronomy, University of Leeds, Leeds, 
LS2 9JT, UK 

Argonne National Laboratory, 9700 S. Cass Avenue, Ar- 
gonne, IL 60439, USA 

* School of Physics, University College Dublin, Belfield, 
Dublin 4, Ireland 

^ School of Physics, National University of Ireland Galway, 
University Road, Galway, Ireland 

Physics Department, California Polytechnic State Univer- 
sity, San Luis Obispo, CA 94307, USA 

Astronomy Department, Adler Planetarium and Astronomy 
Museum, Chicago, IL 60605, USA 

Department of Physics, Purdue University, West Lafayette, 
IN 47907, USA 

Department of Physics, Grinnell College, Grinnell, lA 
50112-1690, USA 

Department of Astronomy and Astrophysics, 525 Davey 
Lab, Pennsylvania State University, University Park, PA 16802, 
USA 

Department of Physics and Astronomy, University of Utah, 
Salt Lake City, UT 84112, USA 

^® School of Physics and Astronomy, University of Minnesota, 
Minneapolis, MN 55455, USA 

Physics Department, McGill University, Montreal, QC H3A 
2T8, Canada 



Department of Physics and Astronomy and the Bartol Re- 
search Institute, University of Delaware, Newark, DE 19716, 
USA 

1" DESY, Platanenallee 6, 15738 Zeuthen, Germany 

Enrico Fermi Institute, University of Chicago, Chicago, IL 
60637, USA 

Department of Physics and Astronomy, University of Iowa, 
Van Allen Hall, Iowa City, lA 52242, USA 

Department of Physics and Astronomy, DePauw University, 
Greencastle, IN 46135-0037, USA 

Department of Physics and Astronomy, Iowa State Univer- 
sity, Ames, lA 50011, USA 

Department of Life and Physical Sciences, Galway-Mayo In- 
stitute of Technology, Dublin Road, Galway, Ireland 

Institut fiir Physik und Astronomic, Universitat Potsdam, 
14476 Potsdam-Golm, Germany 

^® Department of Applied Physics and Instrumentation, Cork 
Institute of Technology, Bishopstown, Cork, Ireland 

^'^ Columbia Astrophysics Laboratory, Columbia University, 
New York, NY 10027, USA 

Author to whom correspondence may be addressed 
^" Deceased 



2 



Acciari et al. 



1. INTRODUCTION 

Gamma-ray bursts (GRBs) have been an active 
area of study since th eir discovery in the late 1960s 
(|Klebesadel et al.l I1973D . Observations of GRBs and 
their afterglows over the last decade are generally consis- 
tent with the relativistic fireball framework (e.g. iPiranI 
Il999f ) . In this theoretical framework, prompt gamma-ray 
emission is produced by internal shocks created by rel- 
ativistic jets with varied Lorentz factors that originate 
from a central engine. The afterglow emission arises 
from external shocks set up when outflowing material 
interacts with the surrounding environment. Within 
this basic fireball framework there have been a num- 
ber of theories proposed that predict VHE photon pro- 
duction. A proposed physical mechanism that produces 
VHE radiation in GRBs is inverse Gompton (IC) scat- 
tering. By this mechanism electrons accelerated by the 
burst's central engine up-scatter relatively soft photons 
fro m an external photon field (external inverse Gomp- 
ton, Bcloborodov''2 005t iWang et al.l l2006^ or from a pho- 
ton field generated by synchrotron emission from the 
electrons themselv es (sy nchrotron self Gompton, SSG, 
iZhang fc Meszarosl[200lt iWang et al.ll2001[ ). 

External shocks may also produce VHE photons. If 
this is the case, measurements of the spectrum above 
10 GeV can directly constrain the medium density as 
well as the equipa rtition fraction of the magnetic field 
(jPe'er fc Waxman|[2005i) in the burst environment. VHE 
emission delayed by ~100 to > 1 0000 s can be produced 
by the external forward sho ck fMeszaro s fc Reei 119941 : 
iDermer et a"I]|2000HFm et a l. 2008). GcV emission from 
electron synchrotron processes in the for ward shock has 
been predicted to be relativel y bright (Zou et al.ll2009D 
and it has been proposed (jKumar fc Barniol DurarJ 
|2009| ) that such emission was detected by t he Large Area 
Teles cope (LAT) on the Fermi satellite (jAtwood et al.l 
120091) in t he bright gamma-ray burst GRB080916G 
()Abdo et~a l. 2009c). In addition to the GeV synchrotron 
component, there may also be SSG p rocesses producin g 
VHE photons in the forward shock (jPanaitescul 120081 ). 
This component is predicted to be less intense than the 
synchrotron component and therefore difficult to detect 
with Fermi-h AT, but the very high energies and rela- 
tively late emission times (up to several hours) make 
these photons prime candidates for detection by ground- 
based, in iaging atmospheric Gherenkov tel escope (lAGT) 
systems ijZou et al.ll2009HXue et al.ll2009[) . 

Yet another possible mechanism for generating delayed 
VHE photons from GRBs is IG scattering of photons 
from X-ray flares. The X-Ray Telescope (XRT) on board 
Swift has made it possible to take detailed X-ray ob- 
servations of fading GRB afterglows on a regular ba- 
sis. In roughly half of these observations. X-ray flare 
activity has been observed that takes place hundreds to 
thousands of seconds a fter the initial gamma-ray signal 
(|Ghincarini et al.l (2007! ). It is predicted that VHE pho- 
tons could arise from the X-ray photons, produced by 
late-time central engine activity, interacting with elec- 
trons accelerated at the forward shock. It is also possible 
that the X-ray flares are produced by the forward shock 
itself and that VHE photons are consequently created 
through the SSG process. Simultaneous observations of 
X-ray and VHE afterglows can distinguish between these 



two possibilities and c an constrain the microphysics in 
the shocks themselves (|Wang et al.ll2006l ). While not ex- 
pected to be a routine event, detection of VHE emis- 
sion from X-ray flares in GRBs by current generation 
lAGTs (VERITAS, MAGIG , HESS) should be possible 
unde r favorable conditions (Fan et a l.l i2008: Galli fc Pir"ol 
120081 ). Recently, the Fermi-L AT detected hard-spectrum 
(F = 1.4) high energy emission associate d with late- 
time X-ray flaring activity in GRB 100728A (jAbdo et al.l 
l20Tl[) . Finally, GRBs have been advanced as a possible 
class of sou rces that generate ultra-high -energy cosmi c 
rays (jWaxmanl 120041 : IMurase et al.l 120081: IDermei1 l200r). 
In hadronic or combined leptonic/hadronic models, VHE 
gamma-rays are produced by the energetic leptons that 
are create d from cas cades initiated by photopion produc- 
tion (Bott cher fc De rmer 1998). 

There have been several attempts to observe VHE pho- 
ton emission from GRBs using ground-based facilities 
but, to date, no conclusive detections have been made. 
A possible detection of VHE gamma rays associated 
with the BATSE-d etected GRB 970417A was reported 
(| Atkins et al.ll2000[ ) by the Milagrito GoUaboration but 
no redshift was determined and no other follow-up ob- 
servations were made. Even though detection of VHE 
afterglow emission with lAGTs is pr edicted to be possi- 
ble, o bservations by both previ ous (iGonnaughton et"al] 
19971) and curren t-generation (jAharonian et al.l 120091 : 
Albert et al.ll2007t) observatories have yielded no signifi- 
cant detections. 

Presented here are the results from GRB observations 
made during an 18 month interval with the Very Ener- 
getic Radiation Imaging Telescope Array System (VERI- 
TAS) between autumn 2007 and spring 2009. The sample 
is limited to well-localized bursts observed with at least 
three of the four VERITAS telescopes. 

2. THE VERITAS ARRAY 

VERITAS is an array of four lAGTs, each 12 m 
in diameter, located 1268 m a.s.l. at the Fred 
Lawrence Whipple Observatory in southern Arizona, 
USA (31°40'30" N, 110°57'07" W). The first telescope 
was completed in the spring of 2005 and the full, four- 
telescope array began routine observations in the autumn 
of 2007. The first telescope was installed at a temporary 
location as a prototype instrument and in the summer 
of 2009 it was moved to a new location in the array 
to make the distance between telescopes more uniform 
and consequently im prove the sensitivity of the system 
(jPerkins et al.l [20091 ). The observations presented here 
were taken with the old array configuration with at least 
three telescopes in the array operational. Each of the 
telescopes is of Davies-Gotton design and is equipped 
with an imaging camera consisting of 499 photomulti- 
plier tubes (PMTs) at the focus, 12 m from the center 
of the reflector. The angular spacing of the PMTs is 
approximately 0°15 resulting in a fleld of view of 3°5 in 
diameter. Each PMT has a Winston cone mounted in 
front of the cathode to reduce the dead space between 
pixels and to increase the light collection efficiency. 

The VERITAS array uses a three-level trigger system 
that greatly reduces the number of background events. 
The first level of the trigger system is at the pixel, i.e. 
PMT, level where the signal from each PMT is fed to 
a programmable constant fraction discriminator with a 
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threshold of 4-5 photoelectrons. The second trigger level, 
the camera/telescope trigger, consists of a pattern trig- 
ger that requires at least three adjacent pixels satisfy- 
ing the first level trigger within a 7 ns coincidence 
window. Finally, an array-level trigger is satisfied if at 
least two of the four telescopes in the array are triggered 
within 100 ns of one another, after correcting for time-of- 
fiight differences. Once the array is triggered, the signals, 
which are continuously digitized for each PMT using 500 
mega-samples per second fiash analog to digital convert- 
ers (FADCs) , are read out and stored to disk. The array 
has an effective area of ~ 10"^ m^ to ~ 10^ m^ and an 
energy resolution of 15 - 20% over the 100 GeV - 30 TeV 
energy range. The single event angular resolution (68% 
containment) is better than 0°.14. A more comprehen- 
sive description of the VERITAS array can be found in 
[Holder et all (120061) . 

3. GAMMA-RAY BURST OBSERVATIONS 

GRB observations take priority over all others in the 
VERITAS observing plan. To facilitate rapid follow-up 
observations of GRBs detected by satellites, VERITAS 
control computers are set to receive notices from the 
GRB Coordinates Network (GCNj!3 over a socket con- 
nection through the TCP/IP protocol. Once the GCN 
notice is parsed by the control computer, an audible 
alarm notifies the observers on duty that a GRB has oc- 
curred. The coordinates of the burst are loaded into the 
telescope tracking software and the observers are noti- 
fied to stop current observations and to begin slewing the 
telescopes to the GRB position, subject to observational 
constraints such as the Moon and horizon. Currently, 
the telescopes are capable of simultaneously slewing at 
a rate of l°s~^ in both elevation and azimuth. Figure [T] 
shows the observation delays for all GRBs with VERI- 
TAS data over a three-year period. The delay between 
the satellite trigger and the beginning of GRB observa- 
tions is usually less than 300 s if the burst is immediately 
observable, and in several cases this delay is less than 100 
s. The dominant contribution to the observation delay 
is the time it takes the telescopes to slew to the source 
position. 

If the GRB is sufficiently well-localized, as is the case 
with the bursts presented here, VERITAS observations 
continue for up to 3 hr after the GRB satellite trig- 
ger, again subject to observing constraints. The tran- 
sition from the prompt to the afterglow phase of a GRB, 
which can occur hundreds to thousands of seconds after 
the initial burst, is often accompanied by X-ray flares 
(jFalcone et al.ll2007D . These flares can be very bright and 
may be associated with extended activity from the GRB 
central engine (Burrows et al. 2005) or be from delayed 
external shocks that could produce a relatively large flux 
of gamma rays in the 100 GeV energy range. For 
GRBs, the VERITAS strategy of rapid follow-up obser- 
vations that continue for several hours allows for good 
temporal coverage of X-ray flare phenomena. Even in 
the absence of flare activity, it is suggested that a sig- 
nificant fiux of high-energy photons from IC processes 
associated with the GRB afterglow may extend to more 
than 10 ks after the beg inning of the GRB prompt emis- 
sion (IGalh fc Pirol[200a) and so an observation window 

http://gcn.gsfc.nasa.gov 



of several hours is warranted. 

During the period beginning 2007 January and end- 
ing 2009 June, VERITAS took follow-up observations 
of 29 GRBs. Nine of these bursts were detected only 
by the Gamma-ray Burst Monitor (GBM) on board the 
Fermi satellite and the errors on the localizations were 
larger than the VERITAS FOV. Analysis of these bursts 
will be presented in a future publication. The remaining 
sample of 20 well-localized bursts (19 detected by the 
Swift satellite and 1 by the INTEGRAL spacecraft) is 
reduced to 16 after applying cuts on the hardware status 
of the array and on the burst elevation (minimum ele- 
vation for GRB observations considered in this analysis 
is 30°). Table [1] lists the general properties of these 16 
bursts. The VERITAS observations of GRBs presented 
here took place during good weather and under dark skies 
or low-moonlight conditions. 

The data were collected in runs with nominal dura- 
tions of 20 minutes with roughly 30 seconds of dead time 
between runs. At the beginning of each run the best 
source localization to arrive via the GCN socket connec- 
tion is used as the target for the duration of that run. 
Twelve of the bursts were observed in "wobble mode" 
in which the source is displaced some angular distance 
away from the center of the camera, allowin g simultane- 
ous o bservation and background estimation (|Berge et al.l 
l2007f) . For the GRB observations presented here, the 
wobble offset was 0%. In the cases of GRB070419A, 
GRB 070521, GRB070612B, and GRB 080604, observa- 
tions were taken in a tracking mode in which the source is 
placed at the center of the camera. Historically, GRB ob- 
servations were taken in tracking mode but wobble mode 
is now the default method of observation with VERITAS 
and all GRB observations are currently taken in this fash- 
ion. The use of the tracking mode does offer a marginal 
increase in "raw" sensitivity over the wobble mode but 
with a significant increase in the uncertainty of the back- 
ground. 

4. DATA ANALYSIS 

The data taken on the 16 GRBs were analyzed us- 
ing t he standard VERITAS analysis suite (jCogan et al.l 
I2008D . The charge in each FADC trace is determined 
by summing the samples over an appropriately-placed 
14 ns-wide integration window. The integrated signal 
from each pixel in the camera results in an image of the 
air shower at the camera plane. The shower image is 
cleaned by eliminating any pixel with a signal of less than 
five standard deviations above its pedestal value, that is, 
a signal less than five times the standard deviation from 
the average FADC measurement when no Cherenkov sig- 
nal is present. Any pixel that registers a signal of at least 
two and a half standard deviations above its pedestal is 
also retained provided it is adjacent to at least one of the 
pixels that exceeds five standard deviations. The cleaned 
images are then parameterized using the Hillas moment 
analysis (pJillas 1985). Before performing a full event re- 
construction, images with less than five pixels surviving 
the image cleaning or with an image centroid more than 
1°43 from the camera center are removed from the anal- 
ysis. A cut on the integrated charge in each image is 
made at ~ 75 (^ 38) photoelectrons for the standard- 
source (soft-source) analysis. For GRBs, the standard- 
source analysis is optimized for a weak Crab-like source 
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TABLE 1 

Details of 16 GRBs observed by VERITAS 
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All information was taken from GCN circulars (http://gcn.gsfc.nasa.gov/gcn3_archive.html). " Duration over which 90% of the emission 
in the 15—350 keV energy band occurs, as measured by the Swift Burst Alert Telescope (BAT). ^ 15—150 keV fluence, as measured by the 
Swift-BAT. T UT time of the GRB trigger determined by the Swift-BAT. ICenko et all ll2007t) . IHattori et"an ||2007| '). HJakobsson et al.l 
II20D7D . 1Prochaska~era l. (2008a). ' iMalesani eTall 1 12008) . IWiersema et al.l ||2008| '). IProchaska et al.l Il2008bl) . ^de Uearte Postigo et"an 
1120091) . HGhornock"era.l., ^OOS l. 
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Fig. 1. — Delay from the start of the burst to the beginning of VERITAS observations for all GRBs with VERITAS data. The open 
symbols correspond to observations that were delayed due to constraints such as the burst occurring during daylight or below the horizon. 
Filled symbols are unconstrained observation delays and are primarily determined by the time it takes the telescopes to slew to the burst. 
The 16 stars correspond to the GRBs discussed in this paper. The shaded regions indicate the annual shutdown of the array due to the 
summer monsoons. 



(3% Crab flux with a spectral index, T = 2.5), while 
the soft-source analysis gives a reduced energy threshold 
and assumes a T = 3.5 spectrum. While the spectral 
characteristics of GRBs are unknown at the highest en- 
ergies, the standard analysis spectral index of 2.5 was se- 
lected based on the average high-energy spectral index, 
/3 observed by BATSE (Kaneko et al. 2006). Since it is 
expected that the extragalactic background light (EBL) 
will significantly soften the intrinsic GRB spectrum, the 
soft source analysis was optimized to the softer assumed 
spectral index of 3.5. It should be noted that although 
the analysis is optimized for a specific spectral index and 
source intensity, this does not preclude the detection of 
sources with characteristics significantly different than 
those assumed. 

At this stage, any event with images in fewer than 
two telescopes is rejected because stereo reconstruction 



is not possible. Furthermore, any event with images in 
only the two telescopes with the smallest separation is 
removed as the proximity of these two telescopes (~ 35 
m) in the old array configuration produced less reliable 
event reconstruction and an increased background rate 
that resulted in decreased sensitivity. After event recon- 
struction, the rejection of background events, which are 
due largely to cosmic rays, is accomplished by comparing 
the length and width parameters of shower images with 
those predicted by Mont e Carlo simulations of ga mma 
ray-initiated air showers ([Krawczvnski et aI]|2006D . Fi- 
nally, a cut on the arrival direction of the gamma ray 
of 61 < 0?13 [0 < 0°.U) for the standard (soft) analysis 
is applied, where 6 is the angular distance in the FOV 
from the reconstructed arrival direction of the event to 
the putative source location. For all bursts presented 
here, the uncertainty in the GRB position (in all cases 
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< 4") is negligible compared to the angular distance cut 
on arrival direction. 

Twelve of the sixteen bursts were observed in "wob- 
ble" mode and the estimation of the background in the 
signal region is made using the reflected region technique 
lAharonian ct al. 2001). In the cases of GRB070419A, 
GRB 070521, GRB070612B, and GRB 080604, the ob- 
servations were made with the GRB positions at the 
center of the fields of view of the telescopes and a re- 
flected region background estimation is not possible. 
For thes e observations th e ring-background estimation 
method (jBerge et al.ll2007D is employed instead. The sig- 
nificance of the gamma-ray excess in the signa l regio n is 
then computed using Equation 17 in lLi fc Mai ()1983D . 

If there is no significant gamma-ray excess detected 
(i.e., the excess in the signal region is less than five 
standard deviations above the background region), the 
99% confidence level upper limit on the number of sig- 
nal photons i s calc ulated using the frequentist method of 
iRolke et al.l (|2005f ). From this number, the correspond- 
ing upper limit on the integral photon flux above the 
threshold energy is computed. The energy threshold is 
defined as the energy at which the product of the detec- 
tor effective area and assumed source spectrum is maxi- 
mized. The effective area, and consequently the thresh- 
old energy, of VERITAS is strongly dependent on the el- 
evation of the source being observed. As the elevation of 
the observation decreases, the column density of the at- 
mosphere increases. This results in a gamma ray of some 
given energy producing a lower Cherenkov photon den- 
sity at ground level, which increases the energy threshold 
of detection. However, because the effective area of the 
instrument is non-zero below the threshold energy de- 
fined in this way, gamma rays in this energy range are 
detectable. For all of the VERITAS data analyzed, a 
secondary analysis was done using an independent soft- 
ware package and the results obtained are compatible 
with those presented here. 

A search for VHE emission is performed over the en- 
tire duration of the VERITAS observations as well as 
over a shorter timescale that optimizes the sensitivity of 
VERITAS to a source with a flux that decays as a power- 
law in time. The Fermi-hAT has detected more than a 
dozen GRBs with emission above 100 MeV. This high- 
energy emission is seen to persist after the flux in the 
GBM band has ceased and shows weak spectral evolu- 
tion with a spectral index between the a and indices of 
the B and function flt to the GBM data (jGhisellini et al.l 
120101) . The temporal behavior of the bright- 
est four Fermi-LA T detected b ursts: GRB 080916C 
(|Abdo et al.ll20093). GRB 090510 (pasauale et al.ll2010l) , 
GRB 090902B~ (|A"bdo et al.l |2009b|), and GRB 090926 A 



(jAckermann et aLll2011t ). show a common ^ t~'^ de 
cay, where 1.2 < A < 1.7 in the observer frame. If 
it is assumed that the temporal and spectral character- 
istics of a GRB detected by the Fermi-LAT extend to 
the VHE energy range, the observed power law tempo- 
ral decay of the high energy emission consequently de- 
fines an optimal duration over which the search for VHE 
emission is maximally sensitive. This optimal duration 
is determined solely by the high-energy temporal power 
law index of the GRB, the delay from the GRB trigger 
time (Ttrig) to the beginning of VERITAS GRB observa- 



tions, and by, to a lesser extent, the observational back- 
grounds. For a VERITAS observation beginning 100 s 
after the GRB Ttrig, the observation window that gives 
maximum sensitivity is ^ 2 - 5 minutes for GRBs sim- 
ilar to the brightest LAT-detected bursts. For bursts 
with unknown high-energy behavior, the determination 
of an optimal time window for VHE observations is not 
straightforward. However, the maximum sensitivity of 
a VHE instrument such as VERITAS to a GRB with a 
power-law decay in time is likely to be on the order of a 
few minutes. 

In the case of GRB 080310, the Swift-XRT detected a 
large X-ray flare beginning ~ 475 s after the beginning 
of the burst as measured by the Swift-BAT. VERITAS 
was on target 342 s after Ttrig for this burst and observed 
throughout the X-ray flare. Figure [2] shows the VERI- 
TAS observ ing window for this bur st relative to the XRT 
fight curve (|Evans et al.l 1 2 0071 [2003) . A search for VHE 
emission is made coincident with the X-ray flare. 

5. RESULTS 

An analysis of VERITAS data associated with the 16 
GRB positions listed in Table [T] shows no signiflcant ex- 
cess of VHE gamma-ray events for any GRB over the 
entire duration of VERITAS observations. Table [5] sum- 
marizes the details and results of the VERITAS GRB 
observations for the sample of GRBs described in Table 
[H The significance distributions for both the standard 
and soft source analyses are shown in Figure |3l The sen- 
sitivity of the VERITAS array, and the small observation 
delays with respect to the GRB Ttrig (half of the burst 
observations had delays of less than 5 minutes) combine 
to give some of the most constraining limits on VHE 
gamma-ray emission from GRB afterglows. 

The VHE photon fluxes from objects at cosmological 
distances are attenuated due to photon absorption by 
the EBL. Of the sixteen bursts for which results are 
presented here, nine had redshifts determined by opti- 
cal followup observations. For these bursts, a limit on 
the intrinsic photon flux of the GRB can be set if one 
assumes a model of the EBL. For all calculations requir- 
in g a model of t he EBL, we use the model described 



Gilmore et al.l ([2009). To determine the factor by 



which the upper limits in Table [2] increase due to ef- 
fects of the EBL, one must calculate the effective atten- 
uation of VHE photons over the VERITAS waveband, 
taking into account the spectral response of the instru- 
ment. For each GRB observation, the effective area of 
VERITAS is multiplied by the assumed intrinsic spec- 
trum of the burst, which we take to be F = 2.5. The 
total flux is then calculated by integrating the intrinsic 
differential flux of the GRB multiplied by the effective 
area of VERITAS, over all energies at which the product 
is non-negligible. This process is repeated, substituting 
an EBL-attenuated burst spectrum for the intrinsic burst 
spectrum. The ratio of the total photon flux obtained us- 
ing the intrinsic burst spectrum to the total photon flux 
obtained using the EBL-attenuated burst spectrum gives 
the attenuation factor for that particular GRB observa- 
tion. For the EBL-corrected upper limits obtained using 
the soft-source analysis, there is an extra correction fac- 
tor to account for the assumed intrinsic burst spectrum 
(r = 2.5) relative to the limits obtained in Table [2] which 
assumes a softer observed spectrum (F = 3.5). The at- 
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Fig. 2.— VERITAS observation window of GRB 080310 superimposed on the Swift-XRT light curve HEvans et ah|[2007l . [2009^ . The inset 
shows the structure of the X-ray flare (dashed lines) and is the time window over which the search for VHE emission was performed. No 
significant excess of VHE gamma rays coincident with the X-ray flare (475 s < t — Tt^ig < 750 s) was found. 




2 3 4 
significance (a) 




2 3 4 
significance (o) 



(a) Standard source analysis 



(b) Soft source analysis 



Fig. 3. — Significance histograms of the 16 GRBs in this sample for both standard and soft source analyses. Included in the figures is 
the normalized Gaussian distribution of mean zero and variance one that the significance histograms should follow if no signal is present. 
The GRB significances are consistent with having been drawn from the aforementioned Gaussian distribution. 



tenuation factors and redshift-corrected upper limits for 
GRBs with known redshift are shown in Table [31 Not sur- 
prisingly, the attenuation depends strongly on both the 
redshift and the energy threshold for a particular obser- 
vation, but under good observing conditions, specifically 
at small zenith angles, reasonable sensitivity out to z ~ 2 
is attainable with VERITAS. 

The search for VHE gamma rays over timescales opti- 
mized for VERITAS sensitivity to a source with ^ ^ 



t^^-^ behavior was performed as described in the pre- 
vious section. Table |3] shows the results of this search. 
No emission associated with any GRB in the sample of 
16 presented in this paper is found. The distributions 
of significances for both the soft and standard optimum 
time analyses are shown in Figure |4l For five of the 
bursts, the maximally sensitive duration of observation 
is greater than the length of time spent observing the 
burst and these bursts are omitted from this analysis. 
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TABLE 2 

VERITAS Observations op Gamma-Ray Bursts 











Standard Source Analysis 


Soft Source Analysis 


o "D 


J- delay (S) 


lobs (mm)^ 


Jilcv. Kangc [ ) ' 


Sth (GeV)* 




Upper Limit 


Sth (GcV)^ 




Upper Limit 


070223 


1.7 X 10** 


74.1 


67-78 


220 


1.3 


9.5 X 10"^^ 


150 


0.8 


2.0 X 10-11 


070419A 


295 


37.7 


32-36 


610 


-0.1 


8.1 X 10"^^ 


420 


-1.0 


1.0 X IQ-ii 


070521 


1118 


75.4 


63-88 


190 


0.1 


4.6 X 10-12 


120 


-0.3 


9.6 X 10-11 


070612B 


201 


131.9 


46-50 


380 


0.6 


2.5 X 10-12 


230 


0.6 


7.1 X 10-12 


071020 


5259 


73.5 


30-43 


570 


1.8 


1.7 X 10"" 


330 


0.5 


2.6 X 10-11 


f\C^ 

080129 


1456 


31.4 


47-50 


370 


1.2 


7.7 X 10^12 


220 


1.4 


1.2 X 10-11 


080310 


342 


198.0 


48-58 


Z/U 


U.Z 


Z.Z X iU 


1/U 


1.8 


(.3 X iU 


080330 


156 


107.8 


64-88 


180 


0.2 


4.0 X 10"" 


120 


-0.7 


6.3 X 10-12 


080409 


6829 


19.0 


31-35 


1300 


0.1 


5.3 X 10"" 


720 


-0.7 


3.8 X 10-11 


080604 


281 


151.8 


33-70 


250 


1.1 


4.7 X 10-12 


160 


0.9 


1.2 X 10-11 


080607 


184 


56.0 


32-46 


400 


1.5 


1.6 X 10-" 


310 


1.1 


2.4 X 10-11 


081024A 


150 


161.2 


55-60 


310 


-1.5 


1.5 X 10-12 


190 


-2.0 


2.2 X 10-12 


090102 


5344 


83.1 


33-48 


400 


-0.1 


8.4 X 10-12 


230 


-0.3 


1.8 X 10-11 


090418A 


261 


30.4 


86-88 


190 


1.0 


1.0 X 10"" 


120 


1.7 


3.0 X 10-11 


090429B 


141 


158.8 


70 88 


180 


1.1 


3.9 X 10"i2 


120 


1.0 


9.6 X 10"i2 


090515 


356 


78.8 


37-57 


340 


0.1 


6.3 X 10-12 


220 


1.2 


2.5 X 10-11 



Upper limits are given at the 99% confidence level in terms of vF,^ at Ethi assuming the spectral indices of 2.5 and 3.5 for the standard 
source and soft source analysis, respectively, in units of erg cm-2 s-l. "Time between the GRB trigger time (Ttrig) and the beginning of 
VERITAS GRB observation. ''Duration of VERITAS observation, t Elevation range of the VERITAS observation. ''The VERITAS energy 
threshold. Statistical significance (standard deviations) of signal counts observed by VERITAS at the GRB position. 



TABLE 3 

ReDSHIFT- CORRECTED VERITAS UPPER LIMITS ON VHE EMISSION PROM NINE Sm/t-DETECTED GRBS 



GRB 


Redshift 


Attenuation 


Standard Source Analysis 


Soft Source Analysis 




(«) 


Factor 


Upper Limit 


Upper Limit 


070419A 


0.97 


1.5 X 10-" 


5.4 X 10-* 


2.8 X 10-® 


070521 


0.553 


0.2 


2.1 X 10"" 


2.9 X 10"" 


071020 


2.145 


1.2 X 10-® 


1.5 X 10-^ 


7.0 X 10-" 


080310 


2.43 


3.1 X lO-"* 


7.0 X 10-^ 


1.4 X 10-® 


080330 


1.51 


0.027 


1.5 X 10-1° 


1.2 X 10-1° 


080604 


1.4 


4.7 X 10-^ 


1.0 X 10-^ 


9.9 X 10-1° 


080607 


3.036 


1.6 X 10-^ 


1.1 X 10"* 


6.8 X lO"-"' 


090102 


1.55 


7.1 X 10-^ 


1.2 X 10-'^ 


8.1 X 10-® 


090418A 


1.608 


0.03 


3.1 X 10-1° 


6.0 X 10-1° 


Upper limit and threshold energy (Eth 


) of each GRB defined 


as in Table 2. The attenuation factor 


is explained in the text. 




significance (a) 

(a) Standard source analysis 




significance (a) 

(b) Soft source analysis 



Fig. 4. Significance histograms obtained from an analysis of the GRBs in the sample over timescales for which VERITAS is maximally 
sensitive to a burst with a t^^-^' power-law afterglow. Both standard and soft source analyses were performed. Included in the figures is 
the normalized Gaussian distribution of mean zero and variance one that the significance histograms should follow if no signal is present. 
The GRB significances are consistent with having been drawn from the aforementioned Gaussian distribution. 
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This occurred when the delay to the beginning of VER- 
ITAS observations was sufficiently long. 

No significant excess of VHE gamma-ray events coin- 
cident with the large X-ray fiare corresponding to the 
interval Ttiig -I- 475 s to Ttiie + 750 s during the afterglow 
of GRB 080310 (see FigureE]) is found. After accounting 
for gamma-ray attenuation by the EBL, the soft source 
analysis returns an integral upper limit of 9.8 x 10~^ 
photon cm^^ s^^ above 310 GeV. Though the flare was 
quite bright in the XRT band, increasing by ~ 3 orders 
of magnitude relative to the underlying afterglow, the 
burst was at a moderate redshift (z — 2.4) so the VHE 
gamma-ray attenuation is presumably significant. 

6. DISCUSSION 

The upper limits on VHE emission presented here pro- 
vide strong constraints on the amount of energy released 
during the early afterglow phase of GRBs. The limits 
themselves, however, are not sufficient to reveal much 
without taking into context the effects of the EBL and 
the intrinsic properties of each GRB. The nine bursts 
with measured redshifts have a mean and median red- 
shift of 1.6 and 1.7, respe ctively. Assuming an EBL 
model (|Gilmore et al.l 12009^. one may convert the up- 
per limits obtained from the VERITAS observations to 
limits on the intrinsic GRB flux as is done in Table [31 
The GRBs without measured redshifts are of less util- 
ity but, as a first approximation, one may assume a 
redshift of z — 2.5, which is the approximate median 
of all of GRB s with known redshift detected by Swift 
(|Gehrels et al.l 2009), to correct for the gamma-ray at- 
tenuation from the EBL. 

After the VERITAS upper limits are corrected for EBL 
effects, we compare the VHE upper limits on the fluence 
above 200 GeV with the fluence of the GRB as mea- 
sured by the Swift- BAT in the 15 - 350 keV energy range 
(jButler et al.ll2007L |2010) that is taken as a proxy for the 
overall intensity of the burst. To account for the different 
delays and durations of the VERITAS observations, we 
calculate tmcd, the time since the beginning of the VER- 
ITAS observations of the GRB at which we expect to 
detect half of the photon signal, assuming a time profile 
of the GRB afterglow of ^ cx t~^-^ that is motivated by 
the high-energy afterglows observed by the Fermi-h AT . 
The ratio of VERITAS upper limit on the fluence above 
200 GeV to the BAT fluence versus imod, is plotted in 
the left panel of Figure [5] for each burst. Since we as- 
sume a known time-dependence of the VHE afterglow, 
we may calculate this ratio for any time period after the 
start of the GRB, which we take to be i — Ttrig > 300 
s. Then for each GRB, we calculate the fractional up- 
per limit on the VHE gamma-ray fluence over the entire 
afterglow (300 < t — Ttrig < oo) relative to the fiuence 
measured by the BAT. A histogram of this quantity is 
plotted in the right panel of Figure O It should be noted 
that if the bursts with unknown redshift are assumed 
to have the mean redshift of the GRBs in our sample 
(z = 1.7) as opposed to mean redshift detected by Swift 
(z = 2.5), then the distribution of bursts with unknown 
redshift moves to the left and more closely follows the 
distribution of known- z bursts. 

These results show that for several bursts the VHE 
component of the GRB afterglow is less than the en- 
ergy released in the Swift-BAT band during the prompt 



phase of the burst. With observation delays often on 
the order of a few hundred seconds, the VERITAS upper 
limits begin to restrict theoretical models in which the 
afterglow from the forward external shock contains an 
SSC component in ad dition to the synchrotron compo- 
nent (IXue et al.ll2"009} ). 

VERITAS observations taken contemporaneously with 
X-ray flares during GRB afterglows are also of interest. 
Over the time period of the flare observed during the af- 
terglow of GRB 080310, the VERITAS upper limits con- 
strain the integral of above 300 GeV to be less than 
9.4 X 10~* erg cm~^ s~^, which is a factor of ~ 12 above 
the peak flux observed by the Swift-XRT in the 0.3 - 10 
keV band. In light of the fact that GRB 080310 was at a 
redshift of nearly 2.5, it is clear that VHE observations 
of a strong X-ray flare from a low redshift GRB could 
challenge some models in which SSC processes produce 
VHE emission simultaneously and with com parable in- 
tensi ty to the X-ray emission during the flare ()Fan et al.l 
120081) and add detail to our understanding of the pro- 
cesses occurring in GRB afterglows. 

7. FUTURE PROSPECTS 

The detection of VHE emission from GRBs in light of 
recent observations by the Fermi and Swift space tele- 
scopes remains a challenging, though not unreasonable, 
prospect. The number of GRBs found by the LAT to 
emit GeV radiation is small, with a detection rate on 
the order of one every few months. Combined with the 
-10% ~ 15% duty cycle of an lACT array such as VERT 
TAS, the probability of simultaneous observation of such 
bursts is not high. On the other hand . >30 GeV emission 
has been detected from both short (jAbdo et al.ll2009a|) 
and long (lAbdo et al.l l2009b f) GRBs and, in the latter 
case, persists well after the prompt phase of the burst. 
Furthermore, these observations indicate that the high- 
energy photon abs orption due to the EBL is not so severe 
(|Abdo et al.ll2010[ ) as to rule out ground-based VHE de- 
tections that i n turn cou ld strongly constrain models of 
GRB physics (jCenko et al. 2011), as well as those of the 
EBL. 

Approximately one of every fifteen GRBs detected by 
the Fermi-GBM is detected by the LAT (provided the 
burst also falls in the LAT FOV). Though they are 
rare, some luminous, LAT-detected GRBs should be de- 
tectable by VERITAS. Taking the spectral and tem- 
poral characteristics of the high energy emission from 
the four brightest Fermi-LAT b ursts: GRB 080916C 
(|Abdo et al.l l2009c^. GRB 90510 (IPasauale et aLllMoh . 
GRB090902B ( Abdo ct 'aLl |2009b|), and GRB 090926 A 
(|Ackermann et al.ii2011ll we estimate the expected flux of 
VHE photons in the energy range of VERITAS as a func- 
tion of time. Figure [S] shows the light curves of three of 
these four bursts that we predict to have been detectable 
by VERITAS. GRB 090510 and GRB090926A produce 
significant signal in the VERITAS band for roughly a 
thousand seconds. GRB 080916C had a redshift of z > 4 
and the VHE emission is extremely suppressed through 
interaction with the EBL. It is observed that even for 
bursts with redshift between 1 and 2, some exceptional 
GRBs may be quite bright in the VERITAS energy band. 
However, in the absence of delayed activity (e.g. that as- 
sociated with X-ray fiares) the power law temporal decay 
of the late-time, high-energy emission necessitates rela- 
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TABLE 4 

A Search for VHE Emission on Timescales Optimized on VERITAS Sensitivity to a Power Law Afterglow Decay ~ t~^-^. 



GRB 


Duratioii(s) 


A^on 


A'ofi 


Standard Source Analysis 

cr" Bth (GeV) Upper Limit 






Soft Source Analysis 
a" Sth (GeV) 


Upper Limit 




2.7 X 10'' 






















07041 Q A 


477 
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14 


0.0 


'yon 


4.U X iU 


2 


42 


— U.9 


A on 


4.6 X 10"^^ 


070521 


1809 


3 


113 


— 1. ( 


1 l\J 


o 1 1 n — 12 
3.1 X iU 


23 


364 


— 0.9 


1 1 A 
110 


1.6 X 10"^^ 


070612B 


325 


3 


21 


0.9 


470 


3.8 X lO^'^'^ 


7 


58 


1.1 


270 


9.3 X lO"'^^ 


071020 


8509 






















080129 


2356 






















080310 


553 


3 


23 


-0.2 


480 


3.2 X 10~" 


13 


55 


1.4 


290 


7.9 X 10"^^ 


080330 


252 





15 


N/A« 


260 


2.4 X 10"" 


6 


43 


-0.2 


170 


1.4 X 10"^° 


080409 


1.1 X 10" 






















080604 


455 


2 


40 


-0.6 


200 


1.5 X 10"" 


9 


128 


-0.3 


140 


3.6 X 10"" 


080607 


298 


4 


16 


1.1 


390 


9.2 X 10"" 


7 


46 


0.3 


250 


1.1 X 10"^° 


081024A 


242 


1 


7 


-0.4 


270 


9.9 X 10~" 


4 


29 





190 


1.9 X 10"^° 


090102 


8647 






















090418A 


422 


3 


16 


0.4 


190 


3.1 X 10"" 


8 


46 


0.4 


120 


6.9 X 10"" 


090429B 


228 


2 


7 


0.8 


200 


9.9 X 10~" 


4 


27 


0.1 


140 


1.5 X 10"^° 


090515 


576 


4 


24 


0.3 


320 


2.7 X 10"" 


11 


72 


0.8 


210 


6.2 X 10"" 



Upper limits defined as in Tabic 2. "Due to the low statis tics, the calc ulation of t he Gaussian significance b y Equation 17 of Li fc Ma 1983 
is not valid. The ratio of Poisson means, as discussed in ;Cousins et al.ii2008i and lZhang fc RamsdenlllQQOI is employed instead, though it 
should be noted that the ratio of Poisson means method is quite conservative in situations with low statistics. '■In the case of zero "on" 
counts, the corresponding Gaussian significance is not defined. 
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Fig. 5. — (a) EBL-corrected VERITAS integral fiuence upper limits above 200 GeV, divided by the fiuence measured by the 5t«i/t-BAT 
in the 15 — 350 keV energy band as a function of tmcd as defined in the text, (b) A histogram of the ratio of the VERITAS integral fiuence 
upper limit above 200 GeV, now integrated over the time period t — Ttrig > 300 s, to the Swift-BAT fiuence. One burst, GRB 080409 is 
omitted in this plot as the VERITAS upper limit on the fiuence of this burst is 11 orders of magnitude above the fiuence measured by the 
BAT. 



lively rapid follow-up observations. VERITAS has made 
several GRB follow-up observations with delays of less 
than 100 s and has a median response time of 328 fE3 
and therefore may be capable of detecting the same high- 
energy component that the Fermi-LAT detects, provided 
it extends to > 100 GeV energies. 

VERITAS continues to take follow-up observations of 
GRBs. In the summer of 2009 one of the telescopes in 
the VERITAS array was moved to a new position that 
resulted in an improvement in sensitivity of ~ 30%. By 
fall 2012, an upgrade of the telescope-level trigger system 
and the replacement of existing PMTs with a more sen- 
sitive PMT will significantly increase the low energy re- 

These numbers are based on all GRBs observed by VERITAS 
from 2007 January through 2009 June, including _Fermi-GBM trig- 
gered observations that are not included in this paper. 



sponse of the instrument. This is particularly important 
for GRB observations as the EBL significantly attenuates 
the high-energy component of sources with appreciable 
redshifts. Additionally, work is ongoing to improve the 
sensitivity of the array with respect to low elevation tar- 
gets, which make up the majority of GRB observations. 
Response times for immediately observable bursts have 
been gradually decreasing and efforts are underway to 
increase the slewing speed of the telescope motors to re- 
duce these times further. Such efforts are increasing the 
GRB science capability of VERITAS and will lead to 
a more thorough characterization of the highest energy 
emission from GRBs. 

8. CONCLUSIONS 

The VERITAS telescope array was used to make 
follow-up observations of 29 satellite-detected GRBs dur- 
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Fig. 6.— Predicted VERITAS light curves for three of the four brightest Fermi-hAT GRBs. The fourth, GRB 080916C h ad a redshift of 
nearly 4.4 and VHE emission is predicted to be too attenuated by the EBL to be detectable by VERITAS. The EBL model of lGilmore et all 
(|2009( ) is used to estimate the attenuation of the VHE 7-rays. The elevation of the burst vfith respect to VERITAS is chosen to be 70° and 
no intrinsic spectral cutoff of the high energy emission is assumed. Each point signifies a detection of at least three standard deviations 
((t) in that time bin. 



ing the period of 2007 January through 2009 June. Due 
to the incorporation of real-time alerts from the GCN 
into the VERITAS pointing and control software, rel- 
atively small observation delays (down to 91 s) were 
achieved. After quality selection cuts, data from 16 of 
the 29 bursts were analyzed and those results are pre- 
sented here. No significant excess of VHE gamma rays 
from any of the bursts is detected and the 99% confidence 
level upper limits on the photon flux are derived. Assum- 
ing a t~^'^ temporal decay of the VHE afterglow, limits 
on the VHE afterglow fiuence relative to the prompt flu- 
ence detected by the Swift-BAT are calculated. For more 
than half of the GRBs with known redshift in our sam- 
ple, the VHE afterglow fiuence is constrained to be less 
than the prompt, low-energy gamma-ray fiuence. 

In the context of recent GRB observations by Fermi- 
LAT, prospects for detection of VHE emission by VER- 
ITAS are good, assuming the most straightforward ex- 
trapolation of the spectral and temporal characteristics 
of the high-energy emission. Contemporaneous early- 



afterglow observations of a GRB by Fermi-LAT and an 
lACT array would provide valuable insights into under- 
standing the physical processes at work in the GRB envi- 
ronment as well as constrain the properties of the EBL. 
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